JOURNAL OF RESEARCH of the National Bureau of Standards— B. Mathematics and Mathematical Physics 

Vol. 64B, No. 2, April-June 1960 

Note on the Solution of Riccati's Differential Equation 

H. Herbert Howe 



(February 4, 1960) 

Three recurrence formulas are developed, giving the solution of a particular case of 
Riccati's equation in power series valid in the neighborhoods of 0, °° , and an arbitrary 
point, respectively. The first two weie programed for computation on the SEAC. 

1. Introduction 



A particular case of Riccati's differential equation 
is important in the theory of radio propagation. 
In the notation often used in that subject, omitting 
subscripts irrelevant in the present paper, the 
equation is 

^_25 2 r + l = (1) 

(It 

where 8 will be considered the independent variable, 
and both variables are complex. The equation 
may be written 

JI = l (2) 

d8 2<5 2 r-l V ; 

Since this is a first-order differential equation, 
there should be one constant of integration. If 
the value of r is arbitrarily chosen for some value 
of 8, the differential equation ought to define r for 
all other values of 8. We may, for example, arbi- 
trarily select t(0) or r (<»), and from it derive r for 
all other points. 

If there exists in the complex 5-plane any point 
P such that as we approach it the value of 8 2 r 
approaches J, then the derivative drjd8 approaches 
infinity, and P is a singular point of r. Since 8 2 t 
involves both the independent and the dependent 
variables, the position of P depends upon the par- 
ticular function r which we are studying; e.g., 
upon the initially selected value r(0). It is not 
immediately obvious that any such P exists, nor 
that there may not be more than one. 

A single-valued function of a complex variable 8, 
with isolated singular points, can be represented 
[1] by power series as follows: In the 5-plane, draw 
concentric circles, one through each singular point, 
and all centered at the origin. Within the smallest 
circle, the function is represented by an ascending 
series in 8 (i.e., there is a constant term, and a series 
of all positive integral powers of 5). This series is 
absolutely convergent and represents the given 
function for all values of 8 within the circle. Outside 
the smallest circle, it diverges. On the circle itself, 
it may or may not converge. This circle is called 
the " circle of convergence" for that series; as a 
special case, its radius may be zero. 



Outside the largest circle, the function is repre- 
sented by a descending series, namely a series con- 
sisting of a constant and negative integral powers 
of 5. It is absolutely convergent outside the circle, 
divergent inside it, and may or may not converge 
on it. 

In the region between any pair of circles, the 
function is represented by a series containing all 
integral powers of <5, both positive and negative, 
and a different series is required for each inter- 
circle region. Such a series is absolutely convergent 
in the region between the two circles to which it 
appertains, may or may not converge on those 
circles, and diverges outside the specified region. 

All of these series together represent the given 
function at all points of the complex plane, except 
perhaps those on the various circles. The function 
exists for all points except the singular points 
themselves; for other points on the various circles, 
one of the series may or may not represent the 
function. 

For the function in which we are interested, we at 
first use the hypothesis that there is one singular 
point, and that therefore there are two series: The 
ascending and descending ones, which have a com- 
mon circle of convergence. A few terms of each series 
are given by Bremmer [2], and a greater number of 
terms are given by Johler, Walters, and Lilley [3]. 

Although one of the two series converges for any 
value of 8 that does not lie on the common circle of 
convergence, the rate of convergence is slower the 
nearer the point is to that circle. For points very 
near the circle, the converging series will not, in 
general, provide a practical method of computation. 
Even within the limits of practical computation, it 
will be impossible to select any fixed number of terms 
at which to stop the computation, if given accuracy 
of the computed result is desired. Rather, the num- 
ber of terms must be increased as 8 gets nearer the 
circle. 

One could, of course, always compute a fixed num- 
ber of terms; but there would then be no assurance as 
to how accurate the result would be, nor as to which 
series should be used, nor even that either one con- 
verged. There is thus need for a formulation of the 
series which does not rely on expressions with a fixed 
number of terms. 
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2. Ascending Series 

In 1955, the present author programed radio appli- 
cations of (1) on the National Bureau of Standards 
Electronic Automatic Computer (SEAC). Some of 
the results have been published [4], but not the for- 
mulas. Instead of using a fixed number of terms of 
the series, recurrence formulas were developed for 
expressing each term as a quadratic function of the 
preceding ones. To derive the formulas for the 
ascending series, we express r as a series in 8: 



r=S M*. 



n=Q 



Kewrite Riccati's equation as follows : 

1 



dr/dd 



-2<5 2 r+l = 0. 



(3) 



(4) 



Insert/into (4) the value of r given in (3), and in one 
term change the summation index from n to m : 



- — 2 s M m+2 +i=o. 



(5) 



n=0 



Now clear of fractions, remembering that two abso- 
lutely convergent series may be multiplied term by 
term: 

1-2 Z) S nb n b m 5 m+n+1 +f] nb n 8 n ~ l =0. (6) 

771 = 71=0 11 = 

In the second term, introduce a new summation 
index k, and eliminate the index m, k being defined by 
m—k—n. To find the summation limits for n and k 
make a plot of m and n; in (6) we are to evaluate the 
term for each point of the following grid: 
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Since k=m-{-n, lines of equal k are the diagonal 
lines. Since we need to sum for all points of the 
grid, extending upward and to the right indefinitely, 
it is seen that we get the same points if we let k run 
from to a, and for each value of k let n run from 
to k. 

l-2Z}i>M*-n^ +1 +i>M w - 1 =0. (7) 

k=0n=0 n=0 

Now change the summation indices again. In 
the second term, replace k by j-2, and in the last 



term replace n by j. Then write out explicitly the 
terms for j = and 1, and collect terms having a 
common summation: 



3=2 L n=0 J 



(8) 



Since this equation is to be satisfied for every 
value of d within the circle of convergence, the 
coefficient of each power of 5 must vanish. This 
imposes no condition on b , which is equal to r(0), 
which we have previously proposed as the constant 
of integration. From the constant term of (8), 



h=-i 



(9) 



and for the higher powers of 8, we see that the 
bracket of (8) must be zero for every value of j^2 r 
i.e., 

bj=(2/j)^nb n bj_ n _ 2 . (10) 

n=0 

In particular, b 2 : =0. 

Now break the right side of (10) into two equal 
sums, each having a coefficient of 1/j; and in one 
of them replace n by a new summation index m 
defined by n=j-2-m: 

1 j-2 I 

^S«M/-»-2+; Z) (j-2— m)6 ; _ 2 _ m & w . (11) 

J n = J m=j—2 

Now replace summation-index m by n. Since the 
summations run over the same values of n, and since 
terms with the same n have also the same 6-product, 
the series can be combined to give: 

&r=^2 b n b^ n . 2 fovj^2. (12) 

J n = 

Equations (9) and (12) show the coefficients to be 
used in (3), after the summation-index of the latter 
has been changed to j. 

3. Recurrence Formulas for the Terms 

For the SEAC computation, it is better to use 
recurrence formulas involving the terms rather than 
the coefficients. Let Bj represent a term of our 
series: i.e., 

Bj=bj8*. (13) 

Substituting into (3), (9), and (12), we get 



r=T J B j ;B 1 =-8;B 2 =0; 
j=o 

Bj= — 8 2 Zj Bnfij-n-'t 

J n=0 



for j>3, 



(14) 



where B o =b =T(0) is the constant of integration, 
and may be chosen arbitrarily. In practice, it is a 
boundary condition, chosen to fit the conditions of 
the problem at hand. 
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4. Descending Series 

Using a similar but somewhat more complicated 
method, we get for the descending series 



(15) 



Ai= {~2Aj\_Si A,,Aj - n 



where 4o=t(qo) is the constant of integration. 

In order to determine which value of A belongs 
to the same function t(8) as does a given value of 
B , so that the results of the two series may be com- 
bined into a single function, other information about 
the problem at hand is used [3, pp. 1 to 4]. 

5. The SEAC Computations 

The SEAC 1 code used (14) and (15) as follows: 
One of them was selected more or less arbitrarily. 
The terms Aj or Bj were computed and saved, since 
each involves those which come before. When 
possible, computations were continued until a term 
was found which was negligibly small. (Real and 
imaginary parts of Aj or Bj each smaller than, say, 
10~ 6 .) This determination began with j=3, since 
B 2 is identically zero. When such a term was 
reached, the computations were discontinued, and 
the A's or B's were summed. 

If any Aj or Bj (starting with A 2 or B z ) was very 
large, that computation was abandoned, and the 
other formula was tried. The 4 same thing was done 
if no negligibly small term was found within the 
space available for storing terms in the memory 
(about 65 terms could be stored). If, by the same 
criteria, the second formula did not converge, the 
computation for that particular value of r(0) was 
abandoned. Some of the gaps in the published 
results [4] come from this cause. (Some occurred 
because of lack of time to complete the lengthy 
computations that had been planned.) This might 
occur either because the point was very near the circle 
of convergence, so that still more terms were needed; 
or because for certain functions there are more than 
one singular point, giving a region in which both 
series actually diverge. 

6. Convergence of the Series 

The criterion used for terminating the computa- 
tions (real and imaginary parts of some term both 
less than 10~ 6 ) is not rigorous, and under some 
circumstances it would be incorrect. This may be 
most easily studied by substituting into (12) so 
that each bj is a polynomial in b =^B , to give the 
equations of [3]. If B happens to be such that any 



bj is zero, then Bj will be zero regardless of 5, even 
though subsequent values of B might be of sub- 
stantial size. For example, two out of every three 
6's are zero if B = 0] or b 7 is zero if 5 + 8/io = 0, etc. 

For the problem under consideration, each value 
of B for which r(<5) was to be derived had an ampli- 
tude of tt/3, and the moduli were all transcendental 
numbers exceeding 1.8 [3, p. 20]. Being transcen- 
dental, they could not make an}- bj exactly zero; 
and it is extremely unlikely that, among the small 
number of values that were used, there Mould be 
one that would make any b 3 negligibly small. The 
expansions of bj through ,7 = 11 as polynomials in 
b [3] show that the coefficients of powers of b tend 
to increase as j increases; and since the modulus of 
b exceeds unity, it seems most probable that the 
values of bj continue to increase as j increases. 
Hence, it seems warranted to conclude that when 
any Bj becomes negligibly small it is probably 
because |<5| is less than unity, and |<5 ; | is small; so 
that succeeding values of \Bj\ probably will be still 
smaller, on account of the factor 8 2 in (14). 

A similar argument seems to apply to the descend- 
ing series. More reliable criteria would be to require 
that each of several consecutive terms of the series 
be negligibly small. 

7. Bridging the Gap 

Although neither series can be used for practical 
computation for points near the circle of convergence 
(or between the circles of convergence, if there is 
more than one singular point), the function itself 
is regular in any region not including a singular 
point; and for a point far from a singular point, 
the derivatives would not be excessive. Hence, if 
we are far from a singular point, it may be entirely 
proper to interpolate across the gap where neither 
series could be used. 

Another method of getting values in this gap is 
by means of analytic continuation. We can use the 
ascending series to give a value of r for some value 
of 8 which we may call q. It is then possible [1, 
p. 196] to represent t as an ascending series in 5 — q, 
which will be valid within the circle centered at q 
and passing through the nearest singular point. 
Such a series, obtained by methods analogous to 
those already described, is 

r=±C,\ C =t(2); C\=-^-~; 
j=o ^C'o2 — 1 

J n=0 

+2g(a- 2 )^g(7 a <7,_ B _ 1 +2 2 S C u O,- u 

J 71 = 71 = 1 

No computations have actually been made with 
this formula. 
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